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Abstract 



i This chapter is intended to relate recent advances in the field of time-frequency signal processing (TFSP) 
to the need for further capacity of wireless communications systems. It first presents, in a brief and heuris- 
tic approach, the fundamentals of TFSP. It then describes the TFSP-based methodologies that are used in 
wireless communications with special emphasis on spread-spectrum techniques and time-frequency array 
processing. Topics discussed include channel modeling and identification, estimation of scattering func- 
tion, interference mitigation, direction of arrival estimation, time-frequency MUSIC, and time-frequency 
source separation. Finally, other emerging applications of TFSP to wireless communications are discussed. 
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25.1 Introduction 

Wireless communications is growing at an explosive rate, stimulated by a host of ™P or ^ 
applications ranging from third-generation mobile telephony, noteless personal commurn, canon s and 
wireless subscriber loops, to rad, o and infrared indoor communica, tons, nomad, c 

tactical military communications. Signal processing has played a key role ,n proving »l— » key 

problems encountered in communications in general, and in wireless commumcatmns m 

important branch of signal processing called time-frequency signal process, ng (TFSP) ha ™™gric> 

the pas, decades |2, 3|. I, provides effective tools fo, analyzing nonstattonary 

content of signals varies in time, as well as for analyzing linear ttme-varymg systems. The purpose 

chapter is to review the methodologies of TFSP applied to wireless communications. 

Fundamental issues in wireless communications include the problems of interference mitigation in 
CDMA (code-division multiple access) or mul.icarrier CDMA (MC-CDMA), and array pr,cesm S for 

M^ng with^l^rneUalling, dtere exist many types (inherent and noninherent) of interference Rising 
degradation in the system performance, hence reducing the system capacity. In addition to the high p y 
resulting from the smart way of providtng multiple access, CDMA has shown advantages over many other 
multiple-access schemes in terms of reducing the effects of interference and mult, pathfadmg I I . C 
has in fact been selected to be the basic building block for third-general, on w.reless common, Catrons 
(wideband CDMA) ( 1 4], and the promising candidate fo, the fourth generation 
CDMA) [ 15] . However, to achieve the best system performance, it is still very crucial to minimize 
of various types of interference in CDMA systems, namely, narrowband interference (NBI) 1 16], multiple 
access interference (MAI) [17], and, for high-data-rate applications, intersymbol interference ^ ^ 

When multiple sensors are available at the receiver side, array processing ec T 
achieve source (mobile) localization or source separation. Source localization is of great importance 
radar/sonar applications bu, also in wireless mobile communications. Mobile locahzatton can add a num- 
ber of important services to the capabilities of cellular systems, including he p for 
emergency services (also known as E-9U problem), location-sensitive billing, fraud protection, and per 

S °Tta probtaTo f source separation or blind source separation (BSS) arises when considering MIMO 
(mlhlLpu! mult. -output, systems where BSS techniques are used to solve the MA, problem and to 

extract the desired information for each user. . • 

' On the other hand, received signals are generally nonstationary due to source signal s nonstationa 

or to the time-varying nature of propagation channels. Indeed, the wireless communication environment 
^ a multipath propagator phenomenon with Doppler effect, where the received signal is not offiy 
affected by active £J S L n01s I b ut a l so by a sum of attenuated, delayed, and Doppler-shifted versions 
of the transmitted signal [5, 6], As a result, the received signals are affected in strength and shape, depen 
m onXent environments (indoor, outdoor, urban, suburban) speed of mobile 

movements and signaling (bandwidth, data rate, modulation, and carrier frequency). The transmitted 
signals undergo serious fading through the propagation channel. Especially in wideband wireless commu 
nk:ations! the^underlying channel er^ibits a random linear time- varying 'TY' m ° St 
commonly assumed to be a wide-sense stationary uncorrelated scattering (WSSUS) process [9-12] 

All these wireless communication problems involve a time-varying context, and therefore 

TFSP techniques should lead to improvements in system performance . , n , f TFSp 

The focus of this chapter is to review such TFSP techniques and clarify the methodologies of TFSP 
that are most relevant for use in wireless communications. The structure of the chapter ^ 
Section 25 2 briefly introduces the basics of TFSP in order to provide a better insight for the review 
the parts that follL Section 25.3 reviews some applications of TFSP techniques in spread-spectrum 
communication systems. Section 25.4 describes array signal processing techniques using time requency 
distributions (TFDs). And finally, Section 25.5 briefly reviews some other issues encountered in wireless 
communications where the TFSP techniques play a central role. 
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25.2 Time-Frequency Signal Processing Tools 

TFSP is a relatively new field comprised of signal processing methods, techniques, and algorithms in which 
the two natural variables time, t, and frequency, /, are used concurrently in contrast with the traditional 
signal processing methods where time and frequency variables are used exclusively and independently. The 
observation of natural phenomena indicates that these two variables, t and /, are usually simultaneously 
present in signals (e.g., natural frequency-modulated (FM) signals such as the song of some birds). Such 
signals are called nonstationary because their spectral characteristics vary with time [2 1 ]. TFSP is designed 
to deal effectively with such signals by allowing their detailed and precise analysis and processing. It also 
enables the design and synthesis of signals and systems with specific time and frequency characteristics, 
suitable to applications such as wireless communications. 

25.2.1 Limitations of Traditional Signal Representations 

The spectrum of a signal (deterministic or random) gives no indication as to how the frequency content of 
a signal changes with time, information that is important when one deals with a large class of nonstationary 
signals such as FM signals. This frequency variation often contains critical information about the signal 
and the process inherent to real-life applications. 1 

The limitation of ‘classical’ spectral representations is better illustrated by the fact that we can find 
totally different signals (related to different physical phenomena), s a (t) and Sf,(f), which yet have the same 
spectra (that is, magnitude spectra) — see Figure 25.2. Representing signals in a way that is useful for precise 
characterization and identification serves as a part of the motivation for devising a more sophisticated and 
practical nonstationary signal analysis tool, which preserves all the information of the signal, and therefore 
discriminates signals in a better way, using one single complete representation instead of attempting to 
interpret magnitude and phase spectra separately. 

25.2.2 Joint Time-Frequency Representations 

25.2.2.1 Finding Hidden Information Using Time-Frequency Representations 

Revealing the time and frequency dependence of a signal, such as a linear FM (LFM) signal, is achieved by 
using a joint time-frequency representation such as the one shown in Figure 25.1. 2 

In this representation, the start and stop time instants are clearly identifiable, as is the time variation of 
the frequency content of the signal described by the linear pattern of peaks. (This information cannot be 
retrieved solely from either the instantaneous power |s(f)| 2 or the spectrum representations |S(/)| 2 . It is 
lost when the Fourier transform is squared in modulus and the phase of the spectrum is thereby discarded. ) 
The spectrum phase contains the actual information about the internal organization of the signal, such as 
details of time instants at which the signal has energy above or below a particular threshold, and the order 
of appearance in time of the different frequencies present in the signal. The difficulty of interpreting and 
analyzing a phase spectrum makes the concept of a joint time-frequency signal representation attractive. 

Example 25.1 

Figure 25.2 illustrates another example of two slightly different signals with the same spectrum that could 
not be properly analyzed without a joint time-frequency representation. Both signals contain three LFMs 
whose start and stop times are different. The differences are not shown easily in the t or f domain, but 



'This information is encoded in the phase spectrum. However, it generally is not used, as it is difficult to interpret 
and analyze a phase spectrum. 

2 In this example the Wigner-Ville distribution (WVD) is used, as it provides the optimal joint time-frequency 
distribution (TFD) for an LFM signal [21]. 
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Fs = 1Hz N = 128 




Frequency (Hz) 




FIGURE 25.1 Time-frequency representation of a linear FM signal: the signal s time domain representation appears 
on the left, and its spectrum on the bottom. 



appear clearly in the t — / representation, which allows precise and simultaneous measurements of actual 
frequencies and their epochs. The B distribution (BD) [18] is used to represent two signals in the t - f 

domain. . 

Figure 25.2 also indicates another significant reason to use joint time-frequency representations ot 

signals: it reveals whether the signal is monocomponent or multicomponent, a fact that cannot be revealed 
by conventional spectral analysis, especially when individual components are also time varying, such as 
the six chirps in the figure. 

25.2.2.2 What Is a Time-Frequency Representation? 

TFSP is a natural extension of both time domain and frequency domain processing that involves repre- 
senting signals in a complete space that can display all the signal information in a more accessible way [2 ] . 
Such a representation is intended to provide a distribution of signal energy, E(t, /), vs. both time and 
frequency simultaneously. For this reason, the representation is commonly called a TFD. 

A concept intimately related to joint time-frequency representation is that of instantaneous frequency 
(IF) and time delay (TD). The IF corresponds to the frequency of a sine wave that locally (at a given time) 
fits the signal under analysis. The TD is a measure of the order of arrival of the frequencies. 

The TFD is expected to visually exhibit in the t — f domain the time-frequency law of each signal 
component, thereby making the estimation of the IF, f(t), and TD, r d (/), easier, and could also provide 
additional information about relative component amplitudes, and the spectral spread of the component 
around the IF (the spread is known as the instantaneous bandwidth, B; (t)). 

The TFD is often expected to satisfy a certain number of properties that are intuitively desirable for a 
practical analysis. Let us denote by p z (t, /) a TFD that is a time-frequency representation of signal z(t). 
p z (t, f) is expected to satisfy the following properties: 

• Pla: The TFD should be real and its integration over the whole time-frequency domain results in 
the total signal energy E z : 




p z (t, f)dtdf = E 



(25.1) 
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Fs = 1Hz N = 256 

Tjrn e-res = 5 B Distribution 
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FIGURE 25.2 Time-frequency representation of two different three-component signals (using B distribution). 



• P lb: It would also be desirable that the energy in a certain region Rin the t — f plane, E Zr , be express- 
ible as in Equation 25.1, but with limits of integration restricted to the boundaries (At, A/) of the 
region R: 

E Zr = f f p z (t, f)dtdf (25.2) 

J At J Af 

which is a portion of signal energy in the band Af and time interval At. 
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P2: The peak of the TFD and the first moment of the time-frequency representation with respect 
to frequency should be equal to the IF of a monocomponent signal: 



m = 



IZo fW 



(25.3) 



P3: For multicomponent signals, the peaks ofthe TFD should exhibit the various IF laws of individual 
components without the nuisance of ghost terms or interferences. 



There are also some other properties, which were earlier seen as strictly needed, but were found later not 
to be, as detailed in Chapter 3 of [3]. For example, early researchers indicated that a TFD should reduce to 
the spectrum and instantaneous power by integrating over one of the variables, so that 




p z (t, f)dt = \S z {f)\ 2 




Pz(t, f)df = |z(t)| 



2 



(25.4) 

(25.5) 



These two conditions are often called marginal conditions. However, it was shown that some high resolution 
TFDs could be generated that did not meet the marginals [3]. 

Another property that was originally seen as desirable is positivity, but it was shown to be incompatible 
with P2 [21], which is more important in practice. 

A number of questions arise from the above: Can we design a TFD that meets the specifications 
listed above? If yes, how can we do it? What are the significant signal characteristics and parameters that 
will impact the construction of a joint time-frequency representation? How do these relate to the TFD? 
How do we obtain them from the time-frequency representation? The answers to these questions are 
important in formulating efficient time-frequency methodologies specifically adapted for applications 
such as wireless communications and will be briefly discussed in the next sections. More details can be 
found in [3, Chapters 2 and 3] and [19]. 



25.2.2.2.1 Physical Interpretation of TFDs 

Most TFDs used for a practical time-frequency signal analysis are not necessarily positive-definite, so they 
do not represent an instantaneous energy spectral density at time t and frequency /. For example, the 
Page distribution describes the rate of change of the spectrum and is defined as the time derivative of the 
running spectrum (spectrum of the signal from — oo to time t), and hence can take both positive and 
negative values: 

o /V 2 

Pz p (t, f) = — / z(u)e~ ,2nfu du 

dt J - 00 

To relate p z (t, f) to the physical quantities used in practical experimentation, we can interpret p z (t, f) 
as a measure of energy flow through the spectral window ( / — A / /2, / + A//2) during the time interval 
(f — At/2, t + At/2). The signal energy localized in this time-frequency domain, (At, A/), is then given by 

rt+At/2 nf+Af/2 

F At.Af = / / p z {t,f)dtdf (25.6) 

Jt- At/2 J f-Af/2 

The larger the window, the more likely £ A1 ; a/ will correspond to a true measure of physical energy. 
The window should be chosen large enough so that At A/ satisfies Heisenberg’s uncertainty relation [2]: 
AtAf > l/(47t). 




25.2.2.2.2 Instantaneous Frequency and Group Delay 

The IF, fi (f ), of a monocomponent signal is a measure of the localization in time of the individual frequency 
components of the signal [2]. 
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The IF, fi(t), of a monocomponent analytic signal 3 z(t) = a{t)e’‘ t>{ ' t) is given by 



m = 



1 dcp(t) 
2i x dt 



(25.7) 



The IF of a monocomponent real signal s (f) = a(t) cos 0(t) is defined as the IF of the analytic signal 
z(f) corresponding to s (t). The expressions of the IF given above do not apply directly to multicomponent 
signals. For such signals, the expression of the IF needs to be applied to its individual components to have 
a meaningful physical interpretation. 

The twin of the IF concept in the frequency domain is called the time delay, r^if). 

The TD of a monocomponent analytical signal z(f) is defined as 



*d( f) = ~ 



1 dd{f) 
In df 



(25.8) 



where 



T{z(t)) = Z(f) = A(f)e^ (25.9) 

T\:\ stands for Fourier transform (FT). Equations 25.7 and 25.8 are similar except for the minus sign, in 
the same way that the FT and inverse FT (IFT) are similar except for a minus sign in the exponent ofthe 
basis function e > 2n l l . 

The TD of a monocomponent real signal s ( t) is defined as the TD of the analytic signal z(f ) corresponding 
tos(r). 

The order of appearance of each time-varying frequency component is the TD. The global order of 
appearance of the frequencies is called the group delay (a mean value of individual TDs). 

The IF and TD describe the internal organization of the signal. Neglecting this information would 
result in a lack of precision of the information characterizing the signal, as illustrated in [3, p. 7, 
Figure 1.1.2]. 



25.2.3 Quadratic Time-Frequency Distributions 

25.2.3.1 Time- Varying Spectrum and the Wigner-Ville Distribution 

To determine why time information appears to be lost when we take the power spectral distribution (PSD) 
and how it can be recovered, let us consider a complex random process Z(f ) with realizations z(f, e ), where 
e represents the ensemble index identifying each realization. To improve clarity, we simply replace z(f, e) 
by z(t). (e is implicit when we say that z(t) is random.) 

The autocorrelation function of z(t ) may be defined in symmetric form as 

R z (t> t) = E {z(t + z/2)z*(t - r/2)} 

= E[K z (t, t)} (25.10) 

where K z (t, r) = z(t + x/2)z*(t r/2) is the signal kernel and £{•} denotes the expected value 

operator. 

The Wiener— Khintchine theorem states that for a stationary signal, the signal power spectral density 
equals the FT of its autocorrelation function. 



3 The analytic signal z(t) associated with the real signal x{t) is defined as 

z(f) = x(t) + jH[x(t)} 

where H{ } represents the Hilbert transform [21], 



